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Abstract: Theoretical predictions for neutrino fluxes indicate that km 3 scale detectors 
are needed to detect certain astrophysical sources. The three Mediterranean experiments, 
ANTARES, NEMO and NESTOR are working together on a design study, KM3NeT, for a 
large deep-sea neutrino telescope. A detector placed in the Mediterranean Sea will survey a 
large part of the Galactic disc, including the Galactic Centre. It will complement the IceCube 
telescope currently under construction at the South Pole. Furthermore, the improved optical 
properties of sea water, compared to Antarctic ice, will allow a better angular resolution and 
hence better background rejection. 

The main work presented in this paper is to evaluate different km 3 scale detector geometries in 
order to optimize the muon neutrino sensitivity between 1 and 100 TeV. For this purpose, we 
have developed a detailed simulation based on the Mathematica software - for the muon track 
production, the light transmission in water, the environmental background and the detector 
response. To compare different geometries, we have mainly used the effective neutrino area 
obtained after the full standard reconstruction chain. 



Introduction 

Most astrophysical source models predict neu- 
trino production rates in addition to those of 
photons and cosmic rays. The goal of neu- 
trino astronomy is to discover and measure 
these high energy neutrinos. Neutrinos offer 
several advantages over traditional astronomi- 
cal messengers. They have a weak interaction 
cross-section, are stable and neutral. So, they 
can propagate through the Universe nearly un- 
affected by the environment and, contrary to 
cosmic-rays, without being deflected by mag- 
netic fields. 

The recent results from gamma astronomy (es- 
pecially the H.E.S.S. telescopes) and theoreti- 
cal predictions for neutrino fluxes indicate that 
at least km 3 scale detectors are needed to de- 



tect certain astrophysical sources. For this 
reason a European consortium including the 
three Mediterranean experiments, ANTARES, 
NEMO and NESTOR is working on a design 
study for a large deep-sea neutrino telescope: 
KM3NeT [1]. This detector will be placed in 
the northern hemisphere (in the Meditcrancan 
Sea) to survey a large part of the Galaxy, in- 
cluding the Galactic Centre. It will so be com- 
plementary with the IceCube telescope cur- 
rently under construction at the South Pole. 

This paper presents NESSY, a new fast and 
detailed semi-analytic simulation based on the 
Mathematica software 1 . Its purpose is to eval- 
uate and optimize different km 3 scale detec- 
tor geometries for neutrino detection in the 

1. Wolfram Research Inc. - www.wolfram.com 
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Mediterranean sea between 1 TeV and 1 PeV. 
The first part of this report describes the 
Monte-Carlo simulations of muon tracks in- 
side the detector and the reconstruction algo- 
rithm used. In the second part, a hexagonal 
detector was used to determine the optimal ef- 
fective area for muon neutrino between 1 TeV 
and 1 PeV as a function of the photomultiplier 
(PMT) density. 

NESSY: a full simulation and anal- 
ysis chain 

When searching for sources, two main param- 
eters are important: the sensitivity is propor- 
tional to the square of the angular resolution, 
while the number of reconstructed events is re- 
flected in the effective area of the detector. To 
determine the best geometry, we need to opti- 
mize these two parameters. A full Monte-Carlo 
simulation and analysis chain has been devel- 
oped with Mathematica. 

This software provides a fast and flexible tool 
to easily scan a large parameter space, in- 
cluding design properties (distance between 
lines and storeys, PMTs configurations) and 
the environmental characteristics (water opti- 
cal properties, bioluminescence background). 

MC generator 

The simulation is organized in three different 
steps. The first step is the generation of muon 
tracks in a volume 3 times the attenuation 
length ( 30 m) of the Chcrcnkov light in water 
around the instrumented volume. The second 
step is the propagation of the light produced 
by the muon and secondary particles induced 
by electromagnetic showers in the sea water. 
The light scattering in the medium is simu- 
lated with an analytic model of diffusion in sea 
water. Photons arriving on a PMT are con- 
verted into hits according to PMT properties 
(transit time spread, quantum efficiency and 
the electronics response). Finally, the last step 
of the MC generation is to add the environ- 
mental background (K40 beta decay and bio- 
luminescence) . 

Several geometries may be considered for the 
future KM3NeT detector including homoge- 
neous or alternative versions such as the ring 
or clustered configuration. Each configuration 



is constituted by vertical lines composed of a 
number of storeys, each containing a number 
of PMTs with given orientation and character- 
istics. 

Selection, reconstruction & analysis 

The first step is a rough selection of the hits. 
Only those in clusters of a minimum of 3 hits in 
adjacent storeys are considered. Furthermore, 
a causality filter relative to the barycentre of 
selected hits is used. The different steps in 
the hits selection procedure results in a purity 
better than 99% of true hits. 

The objective of the reconstruction is to de- 
termine the muon track characteristics - the 
zenith and azimuth angles - using the recorded 
arrival time and the number of hits. Due to the 
non linear dependence of the measurement, an 
iterative method is used to determine the best 
fit of track parameters. The reconstruction al- 
gorithm contains different consecutive fitting 
procedures using the result of the previous one 
as the starting point. The last fit uses a max- 
imum likelihood method for which the proba- 
bility density function was built with the time 
residual histogram obtained from MC data. 

After the reconstruction, the muon effective 
area is calculated from the product of the ge- 
ometric surface and the ratio of the number of 
reconstructed to generated events. To convert 
the effective area for muons to neutrinos, the 
following formula is used. 

A V (6,E V ) =A„{6,E„).R I1 {E I1 ) x N A .a u (E„) 

where A V (Q,E V ) and AJOjE^) are the effec- 
tive area respectively for neutrinos and muons, 
R^(E^) is the effective muon range, Na is the 
Avogadro number, a v {E v ) is the interaction 
cross section, P©(#, E v ) is the absorption prob- 
ability in the Earth and I(E^,E„) is the con- 
version between neutrino and muon energy. 

The figure 1 presents the neutrino effective 
area calculated at 2 different steps: after the 
selection and after the reconstruction. This 
study was made using a hexagonal KM3 de- 
tector composed of 127 lines spaced by 100 m 
(25 storeys each 15 m). The reconstruction effi- 
ciency increases with the energy and its mean 
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Figure 1: Neutrino effective area calculated for 
a hexagonal km 3 detector with different condi- 
tions: without scattering and before the recon- 
struction (blue), without scattering after the 
reconstruction (black) and with the scattering 
after the reconstruction (red). 

value is around 50%. This figure shows also 
the influence of the delayed scattered hits for 
the reconstruction processes on the neutrino ef- 
fective area. Including the scattered hits does 
not change the effective surface. Also, this 
doesn't have a significant effect on the detec- 
tor angular resolution (mean difference of the 
angle between the reconstructed and the true 
muon track) : 0.08 to 0.09° around 5 TeV. So, 
to reduce CPU time, the scattering mode was 
disabled for the following studies. 

Homogeneous and compact geom- 
etry optimisation 

Based on ANTARES technologies, a homoge- 
neous and almost isotropic detector has been 
optimized for neutrino detection capability 
above 1 TeV. The chosen configuration for this 
study is a homogeneous hexagon with 127 lines, 
each with 25 storeys with ANTARES type 
optical modules (3 x 10" PMTs [2]). For 
this work, the environmental parameters of 
the ANTARES site were used. In particular, 
an absorption length of 30 m (integrated over 
the Cherenkov light and the PMT spectral re- 
sponse) and a 40K background of 100 KHz have 
been considered. 

To obtain the best neutrino effective area, we 
have adopted a hierarchical procedure: first 
finding the optimal distance between lines, 
then, using this distance to find the optimal 
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Figure 2: Effect of the distance between lines 
on the neutrino effective area (here, the num- 
ber of lines and of PMTs remains constant). 

distance between storeys then finally finding 
the best PMT configuration in the storey. In 
order to have a fair comparison between differ- 
ent detectors, cuts in the selection procedure 
were ajusted for each tested geometry. For the 
two first steps, the number of lines and PMTs 
in the hexagonal detector remains constant. 

The distance between lines 

Figure 2 shows the neutrino effective area cal- 
culated for different distances between lines (25 
storeys separated by 15 m). The effective area 
is a compromise between the geometrical sur- 
face which grows with the distance between 
lines and the detection efficiency which falls 
with this distance. Indeed, the most compact 
geometry presents the best result below a few 
TeV, but at higher energy, the configuration 
with 100 m between lines has the optimal sur- 
face. The angular resolution and the detection 
efficiency decrease with the energy. The mean 
angular resolution of the detector at around 
5 TeV is 0.06, 0.09, 0.21 and 0.40° respectively 
for the configuration with 50, 100, 200 and 
300 m spacing. 

The distance between storeys 

Figure 3 shows the effect of the distance be- 
tween storeys on the effective area. The op- 
timal distance is 15 m for energy neutrino de- 
tection below 10 TeV and 30 m for higher en- 
ergy. The angular resolution remains constant 
for these two distances. For larger distances, 
the selection procedure becomes inefficient. In 
particular, the time windows to make clusters 
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Figure 3: Effect of the distance between storeys 
on the neutrino effective area (the PMT count 
remains constant in this study, only the height 
of the line is varying) . 



Figure 5: Preliminary neutrino effective area 
optimized with NESSY for the KM3NeT de- 
tector compared to Antares and IceCube [3] 
effective Area. 
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Figure 4: Effect of the PMT count in the storey 
on the neutrino effective area (the lines and 
storeys numbers remain constant). 

of hits become too large for background rejec- 
tion. The neutrino flux of astrophysical sources 
decreases as a power law with the energy, so 
the configuration with 100 m between lines and 
15 m between storeys is found as optimal. 

The number of PMTs per storey 

Figure 4 shows the effective area for the opti- 
mal configuration with 3 and 6 optical mod- 
ules in a storey. Doubling the PMT detection 
surface increases the neutrino effective area at 
low energy only by a factor around 1.6 (1 TeV), 
but does not change it at high energy. As the 
neutrino spectrum of astrophysical sources de- 
creases as a power law with the energy, the 
configuration with 6 optical modules is better, 
but more expensive. 



Summary and Discussion 

Different configurations with a homogenous 
and almost isotropic detector have been tested 
to determine which one gives the optimal re- 
sult for neutrino detection around 1 to 10 TeV. 
A hexagonal geometry with 100 m spaced lines 
and 3 PMTs per storey every 15 m has the best 
effective area and an angular resolution better 
than 0.1°. This result is compared to the Ice- 
Cube effective area [3] and is slightly better 
up to 10 TeV (cf. figure 5). The ANTARES 
performance is also shown in reference. This 
study is the first step to optimize the perfor- 
mance of a km 3 detector. Other parameters 
have to be taken into account such as the envi- 
ronmental properties of the different sites, the 
effect of the energy reconstruction, etc. In the 
KM3NeT consortium, others groups are test- 
ing some non homogenous geometries and have 
obtained comparable results [4]. Moreover its 
performances have to be adjusted to enhance 
the detection of astrophysical sources [5] . 
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